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In this paper we use the conformal properties of the spinor field to show how we can obtain the 
fermion quasi-normal modes for a higher dimensional Schwarzschild black hole. These modes are of 
interest in so called split fermion models, where quarks and leptons are required to exist on different 
branes in order to keep the proton stable. As has been previously shown, for brane localized fields, 
£ — ' the larger the number of dimensions the faster the black hole damping rate. Moreover, we also 

present the analytic forms of the quasi-normal frequencies in both the large angular momentum and 
. the large mode number limits. 
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I. INTRODUCTION 

> 

. With the advent of theories postulating the existence of additional dimensions, there has been much discussion in 
Qs, ' the literature related to the quasi-normal modes (QNMs) of black holes (BHs), for example, in the context of the 
t— I , QNMs of higher dimensional BHs see reference [l[ . By QNMs we refer to the complex frequency modes of oscillation 
which arise from perturbations of the BH, where the real part represents the actual frequency of the oscillation and 
the imaginary part represents the damping due to the emission of gravitational waves. 

Recent investigations of large extra-dimensional scenarios where the hierarchy problem can be shifted into a 
problem of the scale of the extra-dimensions, has led to the somewhat striking prediction that BHs may be observed 
at particle accelerators such as the LHC [4j , for an interesting treatment of mini BHs within an effective field theory 
framework see However, one poignant problem is that in order to suppress a rapid proton decay we need to 
Qh physically split the quarks and leptons. Such models are generically called split fermions models, see, for example, 
reference @. In supersymmetric versions of this idea the localizing scalars and bulk gauge fields will have fermionic 
bulk superpartners. In this respect it is important to consider the properties of bulk fermions. Note that up to now 
only brane localized QNMs have been calculated [8| (where other BH effects have been considered in reference Q). 

That said, the motivations for studying the fermion QNMs from BHs in this paper are two-fold; the first of these 
being from the theoretical point of view, where the lack of any work done in greater than four dimensions 0] with 
Dirac fields is an omission in the literature. Our calculations serve to fill this gap. Secondly, having a complete catalog 
of all QNMs would be a necessary precursor to eventually studying the emission rates of collider produced, or TeV 
scale, BHs. In a forthcoming work we shall present details of BH absorption cross-sections for bulk Schwarzschild 
fermions in c?-dimensions 0]. See also reference [lOj for tense branes in six-dimensions, where bulk fermions on the 
tense brane background defined there have yet to be obtained. 

As such, this paper shall be structured as follows: In the next section we shall discuss how a conformal transformation 
of the metric allows for a convenient separation of the Dirac equation into a time-radial part and a (d — 2)-sphere, in 
d-dimensions. Such a method has already been applied in reference [llT | to the case of low energy s-wave absorption 
cross-sections (which the authors are currently generalizing to higher energies and quantum numbers). However, this 
has not been used yet in the context of BH QNMs. After this we shall present our results for the QNMs, along with 
the analytical forms of the frequencies in both the large angular momentum as well as the large mode number limits. 
Finally, in the last section, we shall make some concluding statements. 
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II. SPINOR RADIAL WAVE EQUATION 



We shall begin our analysis by supposing a background metric which is d-dimensional and spherically symmetric, 
as given by: 



ds z = -f{r)dt 2 + h{r)dr 2 + r 2 dft 2 d _ 2 , 

where dVl 2 d _ 2 denotes the metric for the (d — 2)-dimensional sphere. 
Under a conformal transformation [ill, E3 : 

— _ q2 

$ _> ^ = n-w-u/fy, 

where we shall take VL — 1/r, the metric becomes: 

ds 2 = -^dt 2 + \dr 2 + dfi2_ 2 , 



where tfj — r 



- M-l)/2 



(1) 



(2) 

(3) 
(4) 



(5) 



Since the t — r part and the (d — 2)-sphere part of the metric are completely separated, one can write the Dirac 
equation in the form: 



[(f v t +TVV) <g> i] i>+ W ® (7 a Va) 5d 5 



V = o, 



where (7 s ) 2 = 1. Note that from this point on we shall change our notation by omitting the bars. 
We shall now let Xi be the eigenspinors for the (d — 2)-sphere [T3|, that is: 



(7 a V a ) s _ Xi (±) = ±i (I , , 
where I = 0, 1, 2, . . . . Since the eigenspinors are orthogonal, we can expand ip as: 



d-2 



(±) 



(6) 



(7) 



(8) 



The Dirac equation can thus be written in the form: 



{t'V* +7 r V r +7 5 



±i I 



d-2 



which is just a 2-dimensional Dirac equation with a 7 s interaction. 

To solve this equation we make the explicit choice of the Dirac matrices: 



where the a 1 are the Pauli matrices: 



Also, 



1 

1 



a 2 = 



-i 

1 



a 3 = 



The spin connections are then found to be: 

T t = a 1 



7 5 = (-ia 3 )(a 2 ) = -a 1 



d ( f 



4y/fhJ dr \r 2 



= 0, 



1 

-1 



T r = 0. 



(9) 



(10) 



(11) 



(12) 



(13) 
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From this point on we shall work with the + sign solution, where the — sign case would work in the same way. The 
Dirac equation can then be written explicitly as: 



dt^ \A^Jh) dr\r 2 



r o 9 , i. ... /, d — 2 



9 . '• '< .1+1 



<9r 2\/7 rfr \ r 



ia n + 



d-2 



We shall now determine solutions of the form: 



(.(+) 



A 1,2 ,-iEt ( iG{r) 



/>(+) - 



= 



(14) 



(15) 



where E is the energy. The Dirac equation can then be simplified to: 



VhJ \ IF 



:dG 



d-2\ ( iG 
F 



ia 3 E 



^)(v G 



(16) 



or 



h dr r \ 2 



(17) 
(18) 



It will be convenient to define the tortoise coordinate, r*, and the function W as: 



L± = A. 

h dr dr* 



r \ 2 



(19) 



In which case, our equations can be expressed as: 

d 



dr* 



-W)G = EF, 



dr* 



W]F = -EG. 



(20) 



The equations can then be separated to give: 



d 



dr 2 * 



-—r + V 1 \G = E z G and \- +v 2 )F = E z F, 



d 



dr 2 * 



(21) 



where: 



v l2 = ± d ^ + w\ 

dr* 



(22) 



Since V\ and V2 are supersymmetric to each other, F and G will have the same spectra, both for scattering and 
quasi-normal. Incidentally, for <j>n \ we have these two potentials again. 

Refining our study now to a rf-dimensional Schwarzschild BH, where / becomes: 



f(r) = h-\r) = l-(^-) 



d-3 



and where the horizon is at r = m with: 



d _ 3 frrMr((rf-l)/2) 
h - 7r(d-i)/2(d-2) ' 



(23) 



(24) 



In this case the potential Vi can be expressed as: 



d-3 



At)' 



1 + 4? 



(25) 
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TABLE I: Massless bulk Dirac QNM frequencies (Re(-B) > 0) for a higher dimensional Schwarzschild BH with I > n > 0, 
with d = 5,6, 7, 8, 9 and 10 dimensions. Given the accuracy of the 3rd order WKB method, results are presented up to three 
significant figures only. 



(l,n) Odd d 


d = 5 


d = 7 


d = 9 


1=0 n=0 


725 - 396 i 


1 79 - 809 i 


2 66 - 999 i 


1=1, n=0 


1.32 - 0.384 i 


2.73 - 0.807 i 


3.73 - 1.03 i 


1=1, n=l 


1.15 - 1.22 i 


2.05 - 2.68 i 


2 30 - 3 57 i 


1=2, n=0 


1.88 - 0.384 i 


3.61 - 0.817 i 


4.71 - 1.06 i 


1=2, n=l 


1.75 - 1.18 i 


3.11 - 2.56 i 


3.65 - 3.35 i 


1=2, n=2 


1.56 - 2.03 i 


2.27 - 4.53 i 


1.80 - 6.23 i 


1=3, n=0 


2.43 - 0.384 i 


4.46 - 0.821 i 


5.64 - 1.08 i 


1=3, n=l 


2.33 - 1.17 i 


4.08 - 2.52 i 


4.85 - 3.29 i 


1=3, ii=2 


2.17- 1.99 i 


3.38 - 4.36 i 


3.30 - 5.84 i 


Z=3 ; n =3 


1.96 - 2.84 i 


2.46 - 6.36 i 


1.30 - 8.87 i 


(/, ti) Even d 


d = 6 


d = 8 


d = 10 


1=0, ii=0 


1 28 - 639 i 


2.24 - 0.924 i 


3 05 - 1 05 i 


1=1, n=0 


2.10 - 0.623 i 


3 27 - 936 i 


4.14 - 1.10 i 


1=1, n=l 


1.71 - 2.04 i 


2.24 - 3.18 i 


2.26 - 3.87 i 


1=2, n=0 


2.87- 0.631 i 


4.21 - 0.956 i 


5.14 - 1.14 i 


1=2, n=l 


2.11 - 3.42 i 


3.45 - 3.01 i 


3.75 - 3.59 i 


1=2, n=2 


2.27 - 4.53 i 


2.15 - 5.45 i 


1.26 - 6.95 i 


1=3, n=0 


3.61 - 0.632 i 


5.11 - 0.964 i 


6.08 - 1.16 i 


1=3, n=l 


3.39 - 1.93 i 


4.54 - 2.96 i 


5.06 - 3.54 i 


1=3, n=2 


3.00 - 3.32 i 


3.46 - 5.18 i 


2.95 - 6.340 i 


1=3, n=3 


2.49 - 4.78 i 


2.04 - 7.69 i 


0.327 - 10.0 i 



Next, we shall define, for notational convenience: 



K = I + 



d-2 



A ee r d - 3 {r d - 3 -rjf 3 ), 
with k = i — 1, f , | + 1, . . . . This allows the above potential, V\, to be simplified to: 



(26) 
(27) 



For the case d = 4, this then becomes: 



r 2(d-2) 



kA 1/2 _ r d-3 



d-l 



r d-3 
H 



Vi = (kA 1/2 - r + 3M 



(28) 



(29) 



with A = r(r — 2M), which is just the radial equation for the Dirac equation in the 4-dimensional Schwarzschild BH 
0. 



III. QNMS USING THE IYER AND WILL METHOD 



To evaluate the QNM frequencies we adopt the WKB approximation developed by Iyer and Will [15| , also see refer- 
ences therein. Note that this analytic method has been used extensively in various BH cases [l6[ , where comparisons 
with other numerical results have been found to be accurate up to around 1% for both the real and the imaginary 
parts of the frequencies for low-lying modes with n < I (where n is the mode number and I is the spinor angular 
momentum quantum number). Furthermore, we have also included the n = I modes in our results, but the inclusion 
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of these modes does depend on the number of the dimensions d. The formula for the complex quasi-normal mode 
frequencies E in the WKB approximation, carried to third order beyond the eikonal approximation, is given by [l5| : 



E 2 = [V + (-2Vo) 1/2 A] - i(n + ±)(-2<) 1/3 (l + ft), 



where we denote Vq as the maximum of V\ and 



A 



ft 



(30) 



(31) 



1 



1 



v; 



(4)~ 



1 



2304 \~vf) (67 + 68a2) -288 



(5)' 



1 



(6)' 







, , (19 + 28a 2 ) , 



(5 + 4a 2 ) 



Here 



a = n + — , n = 



0,1,2, - •• , Re(E) > 
-1,-2,-3, ••• , Re(E) < 



and 



v; 



(») _ 



(32) 



(33) 



r»=r.(r max ) 



It is worth mentioning that in the spin-1/2 case it does not seem possible to solve for r max for an arbitrary value 
of d, that is, to find an analytic solution for the roots (unlike the case for fields of other spins [l6(). In the case 
of a bulk spin-0 field and the graviton tensor perturbations on a d-dimensional Schwarzschild background a similar 
analytic expression can be found for r max in d dimensions, see reference (18j |. Thus, we must find this maximum 
numerically using a root finding algorithm. Note that for a given d we can solve for the roots analytically using a 
symbolic computer program, although this is not essential, it does improve the performance of our code. 



IV. LARGE ANGULAR MOMENTUM 



If we now focus on the large angular momentum limit (k — > oo) we can easily extract an analytic expression for the 
QNMs to first order: 



E 2 ^V a - l (n+^)(-2Vo) 1/2 + ..., 



(34) 



where Vq is the maximum of the potential V\, see equation (|28[) . In this limit the potential now takes the form: 



The location of the maximum of the potential is at 



r d-l 



d-l\ d ~ 3 



(35) 



The maximum of the potential in such a limit is then found to be: 

^ ^ 2 2^(d-3) 
k->oo [d-l)T^r 2 H 
In this case we find from the 1st order WKB approximation that: 



E 



2~^d-3 
(d- Vj^^TH 



k — i(n + — ) v 7 d — 3 



(36) 



(37) 



(38) 



This result agrees with the standard result in four dimensions, d = 4, for example see reference jl6| . and is similar 
to the spin-0 result given in reference [l7T |. These limiting values also appear to agree well with the plots made in 
Figure [T] and We have also plotted the QNM n = I = dependence on dimension, d, in Figure [31 
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V. ASYMPTOTIC QUASINORMAL FREQUENCY 



Finally, we can also calculate the quasinormal frequency in the limit of large mode number, n. In this case, as we 
can see from the results in the previous sections, that as n — > oo the imaginary part of E tends to negative infinity, 
ImE — > -co. Hence, we are really looking at the large \E\ limit here. Using the method by Andersson and Howls 
, who have combined the WKB formalism with the monodromy method of Motl and Neitzke [2C| , we make our 
evaluations in this limit. Note that in reference this method has been used to obtain the asymptotic quasinormal 
frequency for the four-dimensional Dirac field. Since we are following the same procedure as in references 0, Ell, 
we shall only show the essential steps, where one can consult references 0, for details. To start we return to 
equation (|2Tj) : 



dr 2 



r, ) G = E 2 G ^ + (E 2 ~ V ± ) G = 0. 



(39) 



Defining a new function, 



A 1 / 2 



(40) 



equation (|39|) can be rewritten as: 



d 2 Z 
dr 2 



R{r)Z = 0, 



(41) 



with 



R(r) 



r 4(d-3) 



E 2 - Vi 



(d-3)(d-2)r d h - 3 (d-3)(rf-l)rf- 3) 



2r 



d-l 



4 r 2(d-2) 



(42) 



The WKB solutions to this equation are: 



±2 



(43) 



where t is a reference point and 



?(r) = R(r)- — 



r 4(d-3) 



E 2 - Vi - —, j 



1_ (d 2 -5d+7)rf l - 3 (d-2) 2 r 2 h (d - 3) 
4 r 2 2r d ~ 1 47- 2 ( rf_2 ) 



(44) 



As \E\ — > oo, the zeros of Q 2 (r), or the turning points, t n , in this WKB approximation, are close to the origin in the 
complex r-plane. In this limit, 



Q 2 (r) 



„d-3 



„d-3 



E 2_(d~2) 2 r 



2^2(d-3) 
h 



4 r 2(d-2) 



(45) 



and the turning points are at: 



Q 2 {r) =0^t 



( H 9 )2 2(rf-3) 
2(d-2) \ a " fJ r h 

AE 2 



(46) 



Asymptotically E is very close to the negative imaginary axis, that is, E ~ \E\e t7T / 2 i as such, the turning points can 
then be represented by: 



(d - 2) 2 r 2 h (d - 3) 
4|£| 2 



l/2(d-2) 



i(2n-l)7r/2(d-2) 



(47) 
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for n = 1, 2, • • • , 2(d— 2). Note that the equation giving the asymptotic quasinormal frequency involves two quantities, 
the first being the line integral from one turning point to the other. Here, we have: 



7 



1 - 



(rf - 2) 2 rf d - 3) 

4£2£2(<i-2) 



-an V2 



1 

2 A 

7T 

2' 



1 



1 



1/2 



(48) 



where we have made the change of variable y = 2E£} d ^/(d — 2)rjf 3 \ The other quantity is the closed contour 
integration around r = r%: 



r=r h 



r d-3 E 



r=r h 



d-3 _ r d-3 



r 



2nir h E 
d-3 ' 



(49) 



The asymptotic quasinormal frequency is then given by the formula, 

e~ 2ir = -(1 + 2cos2 7 ) =► e ^r h E/(d-S) = 1 

'd-3 



=>■ E„ = -i 



2r h 



for large n. In terms of the corresponding Hawking temperature, Tjj = (d— 3)/47rr^, 

E n = —i2t:Thtl. 



(50) 



(51) 



We therefore obtain a vanishing real part for the frequency and find that the spacing of the imaginary part goes 
to 2ttTh regardless of the dimension. Note that these results are in accordance with that of integral spin fields in 
higher-dimensional Schwarzschild spacetimes [2(| [U, [22j ■ 



VI. CONCLUDING REMARKS 



In this paper we have presented new results for the QNMs of a massless Dirac field on a bulk (/-dimensional 
background, see Figure [T] and [2j The results can be compared with the brane- localized results of reference [8J], 
revealing that bulk fermion modes result in much larger damping rates (in both cases larger d results in greater 
damping rates). Some words of caution are necessary, as can be seen for example in the d — 10 result, which plots the 
I — 0, 1, 2 and 3 channels. In the next angular momentum channel, I = 4, the point n = I, crosses the imaginary axis. 
The presence of a branch cut along this axis would force us to choose the positive imaginary value in accordance with 
the constraint equation (|33p . However this does not signify that there are modes which are unstable, but indicates, 
as can be deduced from our asymptotic large n analysis, that the WKB approximation is breaking down. 

Furthermore, if a more detailed analysis such as along the lines of Leaver's approach [23[ confirms the result that 
some of the QNMs have real part equal to zero (algebraically special frequencies), as our WKB results imply for 
ii > I, then such modes would be transparent. However, as mentioned, we cannot say too much based upon the WKB 
approximation in this vicinity. 

Also, given that the BH damping rate increases with dimension d, see Figure[3l we can naively infer that the amount 
of energy available to be radiated as Hawking radiation will increase with d. These issues are left for further study 

0- 

In section IV we investigated the large angular momentum limit and found an asymptotic result that agrees with 
the d = 4 result. In section V we discussed the remaining issue of asymptotic quasinormal frequencies. The expected 
result in c?-dimensions is obtained in such a limit. Note that recently Dirac QNMs have been analyzed under the 
influence of Lorentz violation [24| . It is of interest that the authors introduce a 75 term in their equations similar to 
the one found in our reduced 2-dimensional Dirac equation ©. 

Finally, it may be worth mentioning that split fermion theories often have massive fermions in the bulk, examples of 
this are the higher order modes of the quarks and leptons, which are freed from the localizing potential at high energy. 
Other examples are the bulk Higginos and gauginos in SUSY formulations. In addition to our results, depending on 
the energy scale involved, investigations of these massive Dirac QNMs would also be of interest. 
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FIG. 1: Lines of constant I for massless bulk Dirac QNM frequencies for a Schwarzschild BH in odd d-dimensions. 




FIG. 2: Lines of constant I for massless bulk Dirac QNM frequencies for a Schwarzschild BH in even d-dimensions. 
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FIG. 3: Variation of the QNM frequency with dimension for I = n = 0. Both the frequency of oscillation and the damping rate 
increase with the number of dimensions. 



